Exact and Optimized Perturbative Evaluation of /ic 
for Homogeneous Interacting Bose Gases 



Jean-Loi'c Kneur^'0 and Marcus B. Pinto^^0 



^ Laboratoire Physique Mathematique et Theorique - CNRS - UMR 5825 Universite Montpellier II, Montpellier, France 
Departamento de Fisica, Universidade Federal de Santa Catarina, 88040-900 Florianopolis, SC. Brazil 



In the framework of the 0(N) three-dimensional effective scalar field model for homogeneous dilute 
weakly interacting Bose gases we use the expansion, within the large limit, to evaluate the 
parameter Tc which is directly related to the critical chemical potential fic- This quantity enters the 
order-a^n^'''' coefficient contributing to the critical temperature shift ATc where a represents the 
s-wave scattering length and n represents the density. Compared to the recent precise numerical 
lattice simulation results, our calculation suggests that the large A'' approximation performs rather 
well even for the physical case N — 2. We then calculate the same quantity but using different 
forms of the optimized perturbative (variational) method, showing that these produce excellent 
results both for the finite TV and large- A'' cases. 

PACS numbers: 03.75.Hh, 05.30.Jp, ll.lO.Wx, 12.38.Cy 



An issue that has attracted considerable attention recently and that is associated to the perturbation theory 
breakdown problem is the study of how interactions alter the critical temperature (Tc) of Bose-Einstein condensation 
(BEC). Due to its nonperturbative nature, this is clearly a non-trivial problem. The authors of Ref. P| were able to 
show that the transition temperature for a dilute, homogeneous, three dimensional Bose gas can be expressed at next 
to leading order as 



where Tq is the ideal gas critical temperature, a is the s-wave scattering length and n the number density. A 
similar structure is also discussed in Ref. 0- As far as the numerical coefficients are concerned, the exact value, 
4 = -647rC(l/2)C(3/2)-5/3/3 ~ 19.7518, was obtained using perturbation theory 0]. On the other hand, the other 
two coefficients, ci and C2, are sensitive to the infrared sector of the theory and consequently cannot be obtained 
perturbatively, but they can, through a matching calculation, be expressed in terms of the two nonperturbative 
quantities k and TZ which are, respectively, related to the number density and to the critical chemical potential, 
as discussed in the next section. Some early analytical predictions for Ci included the self-consistent resummation 
schemes (ci ~ 2.90) the 1/iV expansion at leading order (ci ~ 2.33) '3| and at next to leading order (ci ~ 1.71) 
and also the linear 6 expansion (LDE) at second order (ci ~ 3.06) 6]. Today the mostprecise values predict 
ci ~ 1.30 as shown by the lattice (Monte Carlo) simulations which produce ci — 1.29 ± 0.05 [3 and ci = 1.32 ± 0.02 
[a, 1^ whereas very recent analytical studies, up to seven loop, predict ci = 1.27 ± 0.05 (Variational Perturbation 
Theory, and ci = 1.30 (improved LDE, JTj). A clear review of all how those methods apply to the problem can 
be found in Ref. [ij. At the same time, the order-a^n^/'^ non-perturbative coefficient, Cj, has not been considered by 
most authors working on the BEC ATc problem. At this point it is worth recalling that while ci is a function only of 
K, which is related to the number density, C2 is a function of k as well as TZ, which is related to the critical chemical 
potential. This coefficient was first evaluated with lattice simulations which predict C2 — 75.7 ± 0.4 for iV = 2 
and a previous analytical evaluation made use of the standard LDE [TsL [T3| as well as its recently improved version 
PH - The best previous analytical results (in comparison to the lattice ones) were produced with the improved LDE 
which obtains — 73.5 at order-(5^ Very recently, the highest presently available (seven-loop order) perturbative 
contributions and a corresponding new evaluation of TZ within the framework of the above mentioned variational 
perturbation theory, has been performed by Kastening |15|. The purpose of the present work is twofold. First, we 
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use the traditional 1/A'' expansion to evaluate analytically, in the large- iV limit, the above mentioned quantity TZ.. 
This large N result is a useful cross-check of a similar evaluation also presented very recently in Ref. since 
our independent calculation is performed in a quite different way, as discussed in detail below. Using this result 
in conjunction with the results of Ref. ^ for the density number one may finally establish the analytical large-iV 
prediction for at the next-to-leading order. Comparing our large- A'^ results for /ic with the ones produced by 
the lattice simulations we show how the approximation performs very well, already at the first non trivial order, 
even for the homogeneous Bose gas case where N = 2. Second, by considering our 1/iV result as being "exact" we 
perform standard as well as improved (LDE) optimized perturbative evaluations of TZ (or equivalently fic), also at the 
order, so that the convergence and reliability of such a method in the BEC context can be further tested. As 
we recall below, the LDE method essentially introduces a mass term within the effective scalar field model, therefore 
avoiding the above mentioned infrared divergence problems of ordinary perturbation theory. Then, the interest of 
considering LDE in the large- limit is that at the relevant (next-to-leading) order we can define and calculate 
the LDE perturbative series to arbitrary orders for the quantity Tc, and thus check the convergence of the LDE method 
in principle to arbitrary orders. Our optimized perturbation analysis, in the large- A'^ limit, shows how the improved 
versions of this method also perform very well in connection with the ATc problem. 

This work is organized as follows. In the next section we briefly recall the effective model which is considered 
throughout the paper. In Section III we perform the exact analytical evaluation at order of the critical chemical 
potential. Section IV is devoted to the optimized perturbation (LDE) applications both for the large- as well as the 
finite A^ case. We pay special attention to the convergence problem. Our conclusions are presented in Section V. 



II. THE EFFECTIVE SCALAR THREE DIMENSIONAL MODEL 



The studies concerning the equilibrium properties of BEC can be addressed by means of a non-relativistic effective 
theory described by a complex scalar field. In the dilute limit, which is the regime involved in those experiments, only 
two-body interactions are important 1] and one may then consider the following C/(l) invariant finite temperature 
Euclidean action 
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where, in natural units, (3 is the inverse of the temperature, /i is the chemical potential and tua the mass of the atoms. 
At the relevant low temperatures involved in BEC the internal degrees of freedom are unimportant and this can be 
taken as an effective model of hard core spheres with local interactions for which a represents the s-wave scattering 
length. 

The field ip can be decomposed into imaginary-time frequency modes '0j(x, Wj), with discrete bosonic Matsubara 
frequencies Uj = 27rj//3, where j is an integer. Near the transition the chemical potential becomes very small as 
compared to the temperature (|;u| << T) and, since the correlation length tends to infinity, the distances are large 
compared to the thermal wavelength A — \j2'Kfi jmA- Therefore, only the zero modes matter and one may integrate 
out all other modes to obtain an effective three-dimensional scalar action with an 0(2) symmetry, given by |^ 
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where r = —2mA^- and u — ASnamAT. In the large-A^ limit mainly considered in this work, and also in Refs. 
0, IE 0)0) the field cj) in Eq. © is formally considered as having A^ components (0i, i = 1, . . . , A^). In this case, the 
Bose-Einstein condensate effective action Eq. Q is the N = 2 special case of the general 0{N) invariant action. The 
three-dimensional effective theory described by Eq. is super-renormalizable and requires only a mass counterterm 
[A) to eliminate any ultraviolet divergence. Using this action one may use any non perturbative technique to evaluate 
the physical quantities related to the number density 
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and to the critical chemical potential 
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where the subscripts u and in Eq. mean that the density is to be evaluated in the presence and in the absence 
of interactions, respectively, and S(0) is the self-energy with zero external momentum. 

The actual relations between the two nonperturbative coefficients appearing in Eq. Q and these physical quantities 
are given by Q 



and 
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where &2: Eq. Q, is 
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with K2 = -0.13508335373. 

As discussed below, the relation between Tc and S(0) comes from the Hugenholtz-Pines (HP) theorem at the critical 
point. As a matter of fact, when quantum corrections are taken into account, the full propagator for the effective 
three-dimensional theory reads 



G{p) - [p^+r + j:,Up)Y 
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where represents the three-momentum and Srcn(p) represents the renormalized self-energy. At the transition point 
(p^ = 0), the system must have infinite correlation length and one then has 
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This requirement leads to the Hugenholtz-Pines theorem result = — l]i.en(0). Since is at least of order u, it 
would be treated as a vertex in a standard perturbation type of calculation in which G{p) = represents the bare 
propagator. This shows that perturbation theory is clearly inadequate to treat the EEC problem at the transition due 
to the presence of infrared divergences. One must then reccur to nonperturbative methods like the numerical lattice 
Monte Carlo simulations (LS), the analytical 1/N technique or the LDE 17] (see for instance Refs. TS] for different 
types of applications). The problem is highly nontrivial since the Hugenholtz-Pines theorem automatically washes out 
all momentum independent contributions, such as the one-loop tadpole diagrams, which constitute the leading order 
of most approximations. In practice, this means that the first nontrivial contributions start with two-loop momentum 
dependent self-energy terms. However, having reduced the original model, Eq. to the effective three-dimensional 
one, Eq. Q, makes it easier to tackle those contributions since one no longer has the problem of summing over the 
Matsubara's frequencies, which is a hard task when the number of loops increases. 

Here, as stated in the introduction, our aim is to apply the 1/A^ expansion and the LDE to Eq. Q to explicitly 
evaluate the quantity Tc = — 2mAA*c i'm-A being the relevant atom mass), comparing our results with the the lattice 
simulation prediction which, for TV = 2, is r^^ = (0.0019201 ± 0.0000021)m" " 
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III. STANDARD LARGE-iV CALCULATION 



The large- evaluation is carried out, as usual, by assigning a factor of N to any closed loop, and considering 
u ^ 1/N. Then for the massless critical theory given by Eq. (jSJ the first non trivial contributions are at the 
next-to-leading order 1/A^, and come from the tadpole (with bubble insertions) type of graphs (see Figure 1): 
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FIG. 1: Feynman graphs for the tadpole contributions to Tc at order, with the resummed propagator (dotted lines). 




FIG. 2: Feynman graph for the chain contribution to rc at order. 



with the "dressed" (resummed) scalar propagator 



-1 



m-[^ + B{k)] , (13) 



where 



^(^-^)=/(C)^^r(^-^forZ.^3, (14) 

represents the basic one-loop D-dimensional (massless) integral. The evaluation of Eq. can be done exactly 

and most conveniently in dimensional regularization. Actually, up to a trivial factor of — w/6, it is exactly the same 
two-loop integral entering the exact evaluation of (0^) at the 1/A^ order, thus already performed in Ref. [J (see also 
Appendix A of Ref. for details). The result (which is finite in three dimensions) is 

(jV,tad) _N u"^ 

- 6 967r2 • ^^^^ 

The other contributions to rc at the next-to-leading 1/A^ order come from the chain exchange type of graphs (see 
Figure 2): 

,(....0 ^ _s(A^,e.c)(o) = -1 / |A Fik) ^ , (16) 

which can also be calculated exactly in dimensional regularization in a very similar way. The result, however, is 
not finite in three dimensions and has therefore to be regularized and renormalized. Although it turns out to be a 
technically rather standard two-loop calculation in dimensional regularization, we shall give it in some details below 
since, to the best of our knowledge, it has not been done before. 

Due to the non-finiteness of the final result one should use the arbitrary Z3-dimensional result for the propagator 
in Eq. p4l) which reads 
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Taking this result in the (_D-dimensional) integral, Eq. 
one obtains the following integral: 



H16|) . and making a convenient change of variable, k = P^/C-D 
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which we wrote in a form convenient for straightforward dimensional regularization integration. The explicit result is 
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Taking now the limit D ^ 3 and taking into account the standard MS scheme conversion factor: 
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where M is the arbitrary renormalization scale and 7^; ~ 0.51776 the Euler constant one obtains, after renormalization 
by the appropriate counterterm (A) which removes the —Nu'^/{Q x 967r2(3 — D)) divergence, the final result: 
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Numerically, for the standard lattice simulation choice of scale, M = u/3, and the physically relevant number of field 
components, N = 2, this gives r, 
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0.001814945 



Then, adding the two contributions, Eqs. (|15|l and 1)21(1 . one gets 
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which is completely equivalent to the result very recently obtained in Ref. 15] (Eq. (25)), the latter being performed 
in a different independent way. Eq. I|22|) is directly comparable with the available lattice simulation results, since the 
latter were also obtained in the MS scheme. 

As an interesting consistency check, note that the logarithmic scale dependence in Eqs. ((22(1 . ((21(1 is just identical 
to the one already obtained 0] as the coefficient C2 in Eq. (Q: more precisely, putting the appropriate factors 
defined in Eq. 10 to pass from Tc to the relevant coefficients of the expression Q (and also putting back the 
arbitrary TV dependence of the relevant graph instead of its large N value in Eq. 1(21(1 '). we have: (In part) = 
16(7V + 2)7rC(l/2)C(3/2)-5/3/3 ln[A8M/ {N u)]. This is of course expected, since the divergence has a perturbative 
origin and arises solely from the setting sun graph, namely the order in a formal perturbative expansion of expression 
((16(1 (see also Fig. 2). In that sense, the novelty of our result resides more in the non- logarithmic piece of Eq. 1(22(1 . 
as well as the explicitly non-perturbative dependence on the coupling u therein. 



For N = 2 and M = u/3 we obtain numerically the exact large- A*" result r, 



(N) 



0.002166755 7i which, when used 



in conjunction with the value k = —Nu/{967r'^) obtained by Baym, Blaizot and Zinn-Justin |^] gives from Eq. ((7|): 



C2(l/A^ order, A^ = 2) ~ 95.7941 



(23) 



IV. OPTIMIZED PERTURBATION EVALUATION 



In this Section our aim is to test the eventual convergence of the optimized perturbation (LDE method) towards 
the "exact" large- Af result for Vc given by Eq. ((22(1 . This supplements a similar analysis performed in Ref. Tl'l for 
the quantity ci in the large- A^ limit, but for the independent quantity C2, so as to establish the reliability of the finite 
A^ LDE results recently obtained in Ref. jTJ|. According to the usual prescription the LDE can be implemented in 
the BEC case by considering the interpolated action 



-(r-m^) 



4r 



(24) 



6 



which interpolates between a free theory {5 = 0) and the original interacting theory (i5 = 1) jalldillJI- ^'^y physical 
quantity ($) is then evaluated and renormalized exactly at some perturbative order by treating the dummy parameter 
{S) as a small number and m as an arbitrary mass parameter to be fixed later. At the end, one sets 6 = 1. While 
the Sm^ is treated as a new quadratic vertex, the other part is absorbed by the propagator which is now naturally 
infrared regulated by the mass term. Accordingly, the LDE amounts in a first stage to perform an order by order 
perturbative calculation of the relevant Green's functions (here the self-energy contributions 5](0)) but in the massive 
theory as defined by Eq. (|24|l . Thus, any finite order perturbative calculation will produce m-dependent physical 
quantities. To get rid of this unwanted dependence one usually demands that <f> has its optimum value at the point 
where it is less sensitive to m: 

P=0, (25) 
dm 

according to the Principle of Minimal Sensitivity (PMS) to which we shall refer as ordinary PMS (OPMS). 
However, we shall also consider below some specific variations of this method that had been more recently develoDed|lll 
l20j |. which in practice essentially introduce some additional arbitrary variational parameters. As explicitly shown in 
Refs. |2ll l22j| . the whole procedure is compatible with the renormalization program which is implemented exactly as 
in the usual perturbative case. Finally note that, according to the original prescription Q, we have treated r {vc at 
the critical point) as an interaction, since this quantity has a critical value which is at least of order 5. The Feynman 
rules for the interpolated theory, in Euclidean space, are — 5r, 5iv? and —5As, for the quadratic vertices and —5u for 
the quartic vertex while the propagator is given by 

G{p) = [p2+TO2]-i . (26) 



A. Reviewing the perturbative LDE calculation 

To make this paper self contained let us briefly recall how the LDE is employed in the evaluation of Vc by redoing the 
step by step evaluation performed in Rcf. 13] . Here, we shall stay within the large- limit only, while the evaluation 
performed in Ref. |13| was carried out for any value of N . According to the HP theorem, applied in conjunction with 
the LDE Il3j , Sr'^c^ is obtained from the perturbative evaluation of self energy diagrams, S(*)(0). To ordcr-5 one 
has only the plain tadpole contribution and a direct application of the Feynman rules for the interpolated theory and 
dimensional regularization (see Ref. 13] for details) gives, in the large- A^ limit, the finite contribution 

<5r«=-I]W(0) = fa^£ . (27) 

Carrying on to order-(5^, still within the large- A^ limit, one must consider ^3] 

....._.«<o,..4-..^»^!i-.^(«)'^, ,.a, 

where m* = m(l — (5)^^^ must be expanded accordingly. Now, one replaces Svc appearing in the right hand side of the 
above equation, with the value Sri^^ obtained at the previous order so that the right hand side remains of order-(5^. 
Then, one can see how the "double-scoop" contribution (third term on the RHS of Eq. H28|l ) is exactly canceled 
due to the HP condition applied to Tc at first order. It is easy to convince oneself that at any higher order (n) all 
large- A^ contributions of order- {N Su)'^ , with 1 < k < n cancel. Therefore to get non trivial results with the LDE one 
must add the next-to-leading order, 0{N'^ ((5u)*''+^) terms, shown in Figure 1, to Eq. (|27|) . The first comes from the 
"setting-sun" contribution (see left-hand-side of equality in Figure 2) , which has been explicitly evaluated by Braaten 
and Nieto |23 and also in Ref. ^3 giving 
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As already mentioned at the end of section HI, this contribution displays an ultraviolet pole as e = (3 — i?)/2 ^ 
(D — > 3) which is, within dimensional regularization, the only primitive ultraviolet divergence associated with the 



7 



effective super-renormalizable three-dimensional theory, 
mass counterterm coefBcient in the MS scheme 

5Af, = -S- 



The pole produced by this primitive divergence fixes the 
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As usual, this "vertex" must be considered also at higher orders so diagrams whose divergences arise from "setting-sun" 
sub-diagrams maybe rendered finite. One then gets the finite second order result 
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Note that contrary to the original massless case, the massive LDE interpolated action allows for the unique plain 
tadpole contribution, given by Eq. (|27|) . to survive at any order. This 0{N^) term must then appear together with 
all 0{N^^) non trivial terms shown in Figures 1, 2. In fact, as far as the LDE is concerned, the linear term is crucial 

for convergence properties 0|. Finally, the relation between the values of rl^^ evaluated at two different (MS) mass 
scales, Ml and M2, can be readily obtained from Eq. (|31|l and reads 



re(Mi) = r,(M2) -f 



Nu^ /Ml 
18(47r)2 ^ \M^ 
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It is not too difficult to see that this relation will be verified at any order in S. At order-cJ^ the only diagram which 
is scale dependent is the "setting sun" . At a higher order (n > 3) this order-cS^ type of contribution can appear only 
as a subdiagram (I]ssun(p))- At the same order a similar graph appears, but this time Src replaces the setting sun 
insertion. However, the "vertex" 6rc is always replaced by its expansion in S which contains, at order-J^, exactly the 
same scale dependent term as given by the setting sun with zero external momentum (Essun(O)), with a reversed sign 
producing a net combination proportional to Essun(p) — 5]ssun(0) which, in three dimensions, turns out to be scale 

independent |l3j| . This means that, apart from the order-^^ setting sun, all contributions to Sri^'^ are automatically 
scale independent. This perturbative procedure, to order-^** is clearly illustrated in Ref. In the present work 

we are interested in the large- limit only so it is easier to add up the whole series (represented by the diagrams of 
figures 1 and 2) in a closed form and then develop it to the desired order in S as illustrated in the next subsection. 



B. Resumming the large- A'^ series 

As shown in the previous subsection, a main difference between the LDE and the exact large- A^ above calculation 
is that with the former one considers a massive theory, which effectively introduces a mass term in the (dressed) 
propagator. So, for practical reasons the LDE evaluation can be performed as discussed in the previous sub-section 
or alternatively as we do now. Having defined the relevant perturbative series in powers of u for the massive case, we 
can proceed to the usual ordinary PMS (OPMS) or its improved version (IPMS). Accordingly, we apply on the latter 
series the substitutions: 

m —^ m* = m {1 — Sy^"^ ; u —> uS , (33) 

in the standard OPMS-LDE [13, or 

TO^TO(l-(5)^/2[l + (l-a)(5]^/2 ; u-^uS, (34) 

and subsequent generalizations in the IPMS-LDE jllj. The difference between the OPMS and the IPMS is basically 
that within the latter one introduces more arbitrary interpolation parameters at each successive perturbative order 
such that the PMS criterion is generalized by requiring d^/dm = 0, d^^/dm? = 0, etc. It turns out IT] that the 
IPMS essentially reduces a basic problem concerning the OPMS. In fact, as discussed in many BEG applications 
[13I [13 . H^. at high orders the OPMS generates many possible (complex) optimum values. On the other hand, the 
IPMS produces, in general, a unique (real) optimum. In the following we will compare results obtained with both 
methods. 

We shall also consider for completeness another variation of the LDE-PMS introduced in Ref. , that we may call a 
"renormalization group inspired" interpolation. This, more precisely, modifies the simplest LDE substitution formula 
Eq. (|33|) by introducing consistently into the substitution the relevant critical exponent: 



m ^ to(1 - (5)^/"' , 



(35) 
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where lo' = 2 Q./{2 — 77), 77 being the anomalous dimension of the critical propagator ^ l/p'^~^ and ft = f3'{gc), where 
gc is the critical coupling and P{g) the renormalization group beta function (see e.g. Ref. for a review and 

definitions of the critical exponents). The benefit of the latter method in the present case is that the value of the 
critical exponent is known exactly for large N: uj'{N — > 00) = 1. (Incidentally, our more empirical IPMS approach 
introducing more parameters is in practice equivalent to the latter at second order in 5, Eq. (|34|1 . where the only 
additional parameter a plays thus a role similar to the critical exponent to'). In our numerical applications below, we 
shall thus consider those three different forms of the LDE optimized perturbative expansion. 

In order to resum the large- iV LDE type of series one should now evaluate, in place of expression Eq. the 
equivalent contributions coming from the graphs shown in Figure 1 now including a mass term in all propagators, 
including the plain tadpole term discussed in the previous section 
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where the first term on the right hand side represents the plain one-loop tadpole while the series represents the tadpole 
with bubble insertions shown in Figure 1. However, note carefully that at the end the first term of the series will not 
be considered since it does not correspond to any real physical contribution being an artifact of the resummation as 
discussed, for an equivalent series, in Ref. Now, the massive dressed propagator reads 



B{k,m) 
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47rfc 
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(37) 



with k = |k| and similarly for p, q in Euclidean space. Contrary to the corresponding expression in the m = case one 
cannot integrate exactly Eq. H36|) due to the non-trivial dependence in k and m of B(k, m). However, the coefficients 
of the resulting perturbative series at arbitrary orders can be calculated numerically by simple integrals 0, 0| . After 
discarding the spurious term discussed above the result can be written as 
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where the perturbative coefficients are given by 
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with 
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= — arctan — 
z 2 



(40) 



and z — k/m (see e.g. Ref. 0| for more details). Similarly, the contributions of the exchange- type, shown in Figure 
2, are now to be evaluated with the massive propagator B{k,m) in Eq. (|37|l 
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(41) 



Note that the ffist term of this series is nothing but and additional plain tadpole contribution, however with no overall 
A'' factor. This is consistent, since we should in fact consider for our purpose all contributions at the next-to-leading 
order. Explicitly this terms gives a contribution TOM/(127r) to rc, which combines with the first term of Eq. 138|) 
so that the final one-l oop tadpole contribution to Tc is {N + 2)mu/ (247r), which is consistent with the usual result for 
arbitrary iV (see e.g. ■ Next, the second term corresponds to the only divergent and scale dependent contribution 
("setting-sun") discussed above, which may be more conveniently evaluated separately as shown in Refs. [l3ll23l|. All 
other higher order terms are finite and scale independent being easily resummed. Once more, the integral cannot be 
done analytically but arbitrary order coefficients of the corresponding perturbative series are evaluated numerically 
and the final series can be cast into the form: 



18(47r) 
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where 




[A{z)Y . 



(43) 



As expected, the scale dependence (ln(M)) within the perturbative resuh Eq. (|42|) above again has the very same 
coefficient than the corresponding scale dependence of the non-perturbative result, Eq. (|21|l . Note also that our 
selection of graphs takes automatically care of the HP theorem with being appropriately replaced. It is easy to see 
that the order-^^ result obtained in the previous sub-section can be immediately reproduced by our resummed type 
of calculation. The final total 0{N~^) contribution is then obtained by adding Eqs. H38f) and H42f) (with however the 
consistent replacement ^ + 2 in the first one- loop tadpole term of Eq. (|38|l . as discussed above). 



Let us now perform the OPMS-LDE and the IPMS-LDE substitutions discussed above in order to obtain optimized 
results also discussing the convergence properties of each method. Actually, as discussed above in section IV. B, the 
first series given by Eq. H38|l is equivalent, up to a trivial overall factor —u/6, to the series defining the quantity k 
relevant to the AT^/Tc shift in the critical temperature. For the latter series, a good numerical convergence towards 
the exact result has already been obtained in the case of the IPMS On the other hand, the exact analytical 

result for the chain-exchange type of graphs. Fig. 2, is unambiguously defined separately as Eq. |2U- We can thus for 
numerical analysis simplicity discard in a first stage the part of the perturbative contributions Eq. H38I) to concentrate 
on the other non-trivial series Eq. 1)42(1 in order to examine its eventual convergence properties in the framework of 
the OPMS and/or IPMS. Since we shall also consider in a second stage the complete series, it provides in this way two 
numerically independent analysis of the (eventual) convergence properties of the different PMS methods. Accordingly 
the exact result to which the resummation of Eq. H42|l alone should be compared with, is the second term in Eq. H22|l . 
which numerically is equal to 0.00181497 m^. 

We give the results of IPMS and standard OPMS in Table I, second and fourth column respectively. The third 
column indicates the corresponding value of the IPMS leading parameter^ a, cf. Eq. H34|l. Next, we give in the 
last right column the results obtained when considering, following Refs. [iflllOl, the renormalization group inspired 
interpolation method (that we shall dub RGPMS for comparison) , as defined in Eq. H35() . Notice that at second order 
(fc = 2), the IPMS and OPMS result coincide simply because it is not possible to implement an (independent) new 
interpolation parameter, so that the IPMS actually starts at order three. (The empty RGPMS result entry at second 
order is simply due to the fact that for lo' = 1 there is an accidental exact cancellation of the linear in m (tadpole) 
term of Eq. H38|l . such that the RGPMS non-trivial optimization solutions also start at order three). 

One can see that at first rough sight, all different PMS approaches seem to perform very well, all exhibiting a 
numerically convergence behaviour towards the exact result as the order increases. Moreover, all the results are 
already remarkably good at third order. In practice the IPMS method, introducing one more parameter at each 
successive order, becomes quite cumbersome at large orders, so that we somewhat arbitrarily stopped our analysis at 
order twelve, beyond which it became unreasonably time consuming for our calculation (performed with Mathematica 
[2^). When examining in more details the behavior of each different method, we first observe that the (real part of) 
the standard OPMS method results get very close to the exact result at larger and larger orders, letting aside the 
embarassing non-zero but small imaginary parts. However, in Table I we have selected among several (many at large 
orders) real and complex optimal results those closest to the known exact result. Nevertheless, if selecting instead 
only the real solutions we find a slightly less convergent series of approximations, but still quite satisfactory (the real 
solutions are provided in Table I, whenever they exist). Alternatively, the IPMS has the advantage that it gives at 
any order a unique real optimal solution, and seemingly converging somewhat more slowly but regularly towards the 
expected result. As far as the RGPMS method is concerned, some of the optimization results appear at large orders 
to be even better that the two previous ones, being impressively close to the exact one as the perturbative order 
increases. This may be due to the fact that the known large- A^ exact value lo' = 1 enters this interpolation. However, 
the RGPMS suffers a priori from the same drawback than the ordinary PMS, namely that the real optimal solutions 
arc not unique, as explicitly illustrated in the table, and we have also selected the best among many more (mainly 
complex) optimized solutions. It is interesting to remark finally that, as observed similarly in ref. |ll| . as the order 
increases the leading IPMS parameter a seems to approach (though very slowly) a critical value a = 2 (which, upon 



C. 



The large-A'^ case and convergence of the LDE 



^ Similarly to the results in Ref. , for the real IPMS optimum solutions higher order interpolation parameters turn out to be numerically 
much smaller in magnitude than a, so that we omit to give their precise values here. 
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neglecting other higher order IPMS interpolation parameters, would correspond to = 1 in the RGPMS approach). 
Next, as a second stage independent numerical analysis using the same OPMS, IPMS or RGPMS methods, we consider 



TABLE I: IPMS, OPMS and RGPMS results for the large- iV BEG = r^'' at difTerent orders k. r<;^'°'"^' = 0.00181497 ...u^ 



k 




a 




test RGPMS, 2 (^' = 1) 


9 


u.uuzuu 




n nn9nfi 




Q 
O 




1 4"^ 

_L .4:0 


n nni sd^ 


n nni s9q 


4 


n nni s'Hd 


1 4"^ 


n nni roq + 4 i n~^ r n nni r'\7 


n nni si qq 


5 




1.49 


n nnisn9 + i 6 in~^A 


n nni 821 


6 


0.001827 


1.53 


0.001813 ± 9. 10^''/, 0.001837 


0.0018195 


7 


0.001825 


1.57 


0.001817 ±6.6 10"''/ 


0.001816 ± 2. 10"^/, 0.0018189 


8 


0.001824 


1.60 


0.001820 ± 6.7 10"''/, 0.001830 


0.0018155, 0.0018185 


9 


0.001823 


1.63 


0.001820 ± 1.0 10""'/ 


0.0018154 ± 9. 10"*/, 0.0018182 


10 


0.001822 


1.65 


0.001819 ± 1. 10"^/, 0.001827 


0.0018153, 0.0018180 


11 


0.0018213 


1.67 


0.0018199 ±9. 10""/ 


0.0018152 ± 5. 10"*/, 0.0018179 


12 


0.0018208 


1.69 


0.0018199 ±9. 10""/, 0.0018252 


0.00181513, 0.00181776 


20 






0.0018178 ± 1. 10"^/, 0.0018219 


0.00181509, 0.0018173 



the full perturbative series consisting of the sum of expressions defined in Eq. (|38|l and 142f) . respectively. The 

results of the optimization are shown in Table II. The numerical behaviour appears to be similar to the one in Table I, 
namely with all different PMS methods appearing to converge quite well. However, we stress once more the uniqueness 
of the real IPMS optimized solutions at each orders, exhibiting thus a relatively slower but completelly unambiguous 
apparent convergence for this method. Finally, if neglecting the imaginary parts of the optimized solutions (which are 
however very small in this case), the RGPMS approach appears again to give the best results at large orders, though 
real solutions do not always exist at arbitrary orders, just like in the OPMS case. 



TABLE IL IPMS, OPMS and RGPMS results for the large-iV BEG rT at different orders k. ri'^^ = 0.00216675 ...u'^ 



k 


IPMS 1 2 


a 


best PMS 1 2 


best RGPMS 1 2 r 1 ■,\ 
rc /u (UJ = 1) 


2 


0.002063 




0.002063 




3 


0.001973 


1.14 


0.001960 ± 7.9 10""/ 


0.002060 


4 


0.002041 


1.70 


0.001991 


0.002101 ± 2.1 10"^/ 


5 


0.002065 


1.73 


0.002020, 0.001995 


0.0021254 


6 


0.002081 


1.75 


0.002033 ± 3.4 lO"''/, 0.001999 


0.002135 ± 7.5 10""/ 


7 


0.002093 


1.76 


0.002046, 0.002001 


0.0021422 


8 


0.002101 


1.78 


0.002055 ± 1.6 10""/, 0.002002 


0.0021463 ± 3.7 10""/ 


9 


0.002108 


1.79 


0.002067 ± 1.7 10""*/, 0.002061 


0.0021495 


10 


0.002114 


1.80 


0.002076 ± 1.7 10""/, 0.002061 


0.0021517 ±2.2 10""/ 


11 


0.002118 


1.81 


0.002083 ± 1.7 10""/, 0.002071 


0.0021535 


12 


0.002122 


1.82 


0.002089 ± 1.8 10""/, 0.002006 


0.0021549 ± 1.5 10""/ 


20 






0.002121 ± 1.3 10""/, 0.002008 


0.0021605 ± 7 10"^/ 



D. The finite iV = 2 case 



Here, for completeness, we extend some of the results of Ref. for rc at arbitrary A'' by incorporating the very 
recently available improved six-loop and seven-loop relevant perturbative contributions performed in Refs. |l(lll5il2til |. 
One obtains to order-u^ and for = 2 
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.(5) 



6tt 





1 




.^^ ( "1 ) 








+ 2- 





j=3 



m 



(44) 



where in our normalization, A2 = 1.40724 x 10"^^ ^ 8.50891 x 10"^ A4 = 3.46941 x 10"^, A5 = 2.23296 x 10"^ 
Aq = 1.81417 X 10-^ and Aj = 1.71627 x 10"^. Performing the IPMS-LDE and OPMS-LDE substitutions, aheady 
discussed, one obtains the optimized LDE results shown in Table III below ^ for M = u/S. As one can see, the 
results of our alternative IPMS method show a remarkable agreement with the numerical lattice result. Concerning 
the results of the standard OPMS method, as shown in the right column of Table III, they approach closely the 
lattice result at third order while the fourth order result departs from the LS value. Then, the fifth and higher orders 
prediction approache the LS result again, considering again here the real parts of complex optimal solutions (although 
the imaginary parts are very small in this case). 



TABLE III: IPMS versus standard OPMS results for the physical N = 2 BEG rc, and correspondingly c" at different orders k. 
r^^ = (0.0019201 ± 0.0000021) 



k 




c'i{IPMS) 


a 




2 


0.00315 


101.2 




0.00315 


3 


0.001828 


69.85 


1.77 


0.00215 ± 7.8 lO'^I 


4 


0.0018988 


78.31 


2.02 


0.00247 


5 


0.0018992 


74.44 


2.116 


0.00226 ±3.3 10~"*/ 


6 


0.0019034 


75.51 


2.183 


0.00205 ± 4.5 10"*/ 


7 


0.0019026 




2.237 


0.001887 ±2.0 lO""/ 



V. CONCLUSIONS 



Considering an effective three-dimensional model for homogeneous dilute weakly interacting Bose gases we have 
evaluated the physical quantity which is directly related to the critical chemical potential. Our evaluations made use 
of the standard 1/7V approximation as well as different variations of the optimized perturbation (linear 6 expansion) 
method. As discussed in the text, this quantity is important in the determination of a numerical coefficient (C2) 
appearing at order-a^rt^/"^ in the expansion of ATc- Surprisingly, the large- iV result rv ~ 0.00216672 compares well 
with the most recent lattice simulation result, r^^ = (0.0019201 ± 0.0000021) even for the case N = 2. As a 
matter of fact the large- iV approximation seems to work better for the determination of rc than for the {(f)'^) prediction 
performed in Ref. Q. Then, having established the large- result we have used the linear 6 expansion, in the large- A^ 
limit, to examine its convergence properties, in a way similar as it was done in a previous work 14] for the quantity 
{(f)'^). Our analysis shows that the standard optimization procedure produces converging results which can be better 
controlled if one uses a recently improved version of the method |llj |. In this way one is able to obtain, at each order 
in (5, only one real optimized solution, contrary to what happens if one uses the standard optimization procedure. We 
also obtained that another variation JOJ of the LDE, introducing consistently the critical exponent lo' = 1 performs 
very well in the large A^ limit. Apart from determining the large- A^ result our work shows the reliability of the LDE 
finite A^ results for Vc obtained in Ref. plf. Finally, considering the large- AT |3| and IPMS-LDE [uj results for n, 
one can establish the analytical values for Cj. For N = 2 they are C2 — 95.79 (large-A^ approximation). Also, our 
numerical results in Table III for C2 from the IPMS-LDE at the latest available orders {5^, S^) compare very well with 
the lattice simulation result C2 = 75.7 ± 0.4 

Note added: After completing a first version of this paper (eprint archive: cond-mat / 0408538| vl) we became aware 



^ The slightly different results at orders 3-5 of the IPMS method for with respect to Table VIII of Ref. are simply due to our 
updated use of the most recent analytical and numerical improvements in the determination of the perturbative series coefficients given 
in Ref. 15]. Note that Cj cannot be determined at order 5^ since this requires the use of some eight- loop contributions to ((/>^) wrhich 
are not available at present. 
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of the very recent work by Kastening '15'|, where an analytical evaluation of rc in the large N limit is also performed 
as well as the relevant perturbative contributions up to seven loop for arbitrary N (partly using some previously 
unpublished seven- loop material by B. Nickel), with a corresponding new evaluation of rc in the framework of the 
"variational perturbation theory" (VPT). Our large N exact result for rc, Eq. H22|l . maybe be considered as an 
independent cross-check of his result, since ours is performed with a different method. Concerning his numerical VPT 
results for = 2, at seven-loop order they lie in the range'^: rc/u^ ^ 0.0019028 — 0.0019401, thus appearing very 
consistent with our results in Table III. 
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